arXiv: 1509.00301 v2 [gr-qc] 22 Oct 2015 


Randers geometry as MOND/dark 
matter 


Qasem Exirifard 

Institute for Research in Fundamental Sciences (IPM), 

Tehran, Iran 

E-mail: exir@theory.ipm.ac.ir 

Abstract. We consider a deviation of the physical length from the Riemann geometry 
toward the Randers’. We construct a consistent second-order relativistic theory of gravity that 
dynamically reduces to the Einstein-Hilbert theory for the strong and Newtonian gravity while 
its weak gravitational regime reproduces MOND and the gravitational lensing attributed to 
the dark matter halo. It also naturally accommodates the observed value of the cosmological 
constant. We show that it predicts a few percent deviation for the post Newtonian parameter 
7 in a part of the regime that interpolates the Newtonian regime to the MOND regime. 
The deviation is consistent with the reported observations but can possibly be detected by 
fine-tuned refinements of the current data or specified future observations. 
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1 Introduction 

MOND was proposed in 1983 by Mordehai Milgrom [1], It was, then, altered to the AQUAL 
theory [2], a non-covariant modified theory of gravity capable of reproducing MOND. TeVeS, 
as a relativistic MOND theory, was proposed by Bekenetsein in 2004 [3]. TeVeS theory, despite 
its successes, does not reduce to the Einstein-Hilbert gravity for strong gravity. The purpose 
of this research is to introduce a relativistic MOND theory that reduces to the Einstein- 
Hilbert gravity in the strong regimes of gravity. In so doing we introduce a deviation of the 
Riemann geometry toward Finsler/Randers geometry as the cause of dark paradigms. 

We adapt the Randers geometry as the physical geometry in the section 2. The dynami¬ 
cal degrees of the freedom in the adapted geometry are scalars, gauge fields and a metric. We 
introduce two gauge fields and two scalars. This introduction requires four coupling constants, 
a functional and a length scale. 

In the section 3 we fix the dynamics such that the strong gravity limit of the theory co¬ 
incides to that of the Einstein-Hilbert gravity. This reduces the number of coupling constants 
to two. This causes one of the gauge fields and one of the scalars to have the kinetic energy of 
ghost. We prove that the classical solutions coincides to those of the Einstein-Hilbert theory. 
We also prove that the ghosts are decoupled from the physical spectrum by performing an 
exact integration over the gauge fields and scalars in the Feynman path integral. Notice that 
this does not constrain the dynamics of the fields in the weak regime of gravity/fields. We 
shall choose the dynamics such that the theory reproduces MOND. 
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The section 4 reports the static solution in the Newtonian and MOND regimes of the 
theory. The sub-section 4.1 calculates the light trajectory around a static solution. We require 
the 7 -the post-Newtonian parameter of theory- to be 1 in the deep MONDian regime. This 
requirement reduces the coupling constants of the theory to one constant and maps the length 
parameter of the theory to the critical acceleration in the MOND. 

We consider a family of the interpolating functionals in the section 5. We ignore the 
negligible contribution of the gauge fields and scalars to the metric in the MOND regime. 
We then find the exact static spherical solution of the theory. We calculate the 7 parameter 
around spherical static solutions in whole of the space including the regime that interpolates 
the Newtonian regime to the MOND regime. We report that requiring the 7 parameter in 
the interpolating regime to be consistent with ref. [5] sets a strong upper bound on the 
coupling constant of the theory as depicted in the fig. 3. We observe that these bounds are 
in agreement with the value of 7 measured at the surface of the Sun as depicted in fig. 4. 

We study the Lunar system constraints on the theory in the section 6 . We report that 
Laser Lunar Ranging measurements [ 6 ] combined with the LAGEOS data [7] provide a lower 
bound on the coupling constant of the theory, as depicted in fig 5a and 5b for various values 
of p that encodes the precision of the used perturbation theory around the Earth. 

We study the consistency of the spherical static solution in the section 7. We see that 
as we move toward the asymptotic infinity, the neglected contribution of the gauge-fields and 
scalars to the energy-momentum tensor accumulates in the deep-MOND regime. For any 
given mass, we calculate the length scale wherein the energy-momentum tensor of the gauge 
and scalar fields must be considered, the result is expressed in (7.18). We report that for any 
given astronomical mass, this length scale is far larger than the size of the observed Universe. 
This is one more illustration that the ghosts do not affect the classical solutions. We do not 
study how to quantise gravity that here is composed of a metric and two scalars and two 
gauge fields. We suffice to the study of the classical solutions (1R solutions) and we show that 
the theory is consistent in its 1R limit. Since the strong held limit of the theory coincides 
to the Einstein-Hilbert theory we expect its UV limit to be the same as the UV limit of the 
Einstein-Hilbert theory. 

There is a freedom in identifying the constant of the integration in the functional of 
the theory. We use this freedom and naturally accommodate the observed value of the cos¬ 
mological constant in the theory in the section of 8 with a tuning at the order one. This 
accommodation identifies the only free coupling constant of the theory. The identified theory 
predicts that the 7 parameter in the interpolating regime deviates from 1 by few percents 
only in a section of the interpolating regime. This prediction can be attested by refining the 
current data or specified future observations. 

We study the cosmology of the theory in the section 9. We slightly change the coupling 
of the gauge fields to the matter current in order to have a homogeneous isotropic time- 
dependent solution. This change leaves intact the static solutions of the previous sections. 
We report that the isotropic homogeneous solution of the theory coincides to that of the 
Einstein-Hilbert theory with a cosmological constant. So the late time cosmology of the 
theory is in accord with observation. 

We discuss on the results and concludes in the section 10. 
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2 Adapting the Randers geometry 

We assume that the orbit of a particle of mass m is derived from the variation of 


,S[x(t)] = — m 


dr(J e 2( t , g fll/ x' J -x 1 ' + A^) 


( 2 . 1 ) 


where 4> is an scalar, A^ is a gauge held and g^ w represents the metric. This is the definition 
of length in the Randers geometry [4]. Notice that the hat metric in our notation reads: 

7 /^ = ( 1 ,- 1 ,- 1 ,- 1 ). ( 2 . 2 ) 

Now let a positive auxiliary held r/(r) be defined on the world-line of the particle and consider 

S[x(t), i)(t)] = -m J dr (^-^e 2 ^g^ u x fi x u + + A^xA , (2.3) 

the variation of which with respect to r/(r) yields 

SS - 

— = o -> t](t) = S’g^x^iy , (2.4) 

inserting which into (2.3) yields (2.1). So (2.3) and (2.1) are equivalent. (2.3) is invariant 
under the reparameterization of the world-line: 

r —7 f = f (r) (2-5) 

(It 

v(t) -> V(r) = ( 2 . 6 ) 

This reparameterization allows to set 


V(r) = 1, (2.7) 

which, hereafter, is referred to as the standard parameterization. The standard parametriza- 
tion simplifies (2.3) to 

S[x(t )] = -m J dr x v + A^x^ , (2.8) 

where 

e^g^x^x" = 1, (2.9) 

is understood as the gauge-fixing constraint imposed on its equations of motion. Notice that 
(2.8) is equivalent to (2.1) but it is simpler. We will make extensive use of this simpler 
Lagrangian througout this work. 1 

lr To string theorists, the difference between (2.8) and (2.1) is similar to the difference between Numbo- 
Goto action and Polyakov action. We could have worked with (2.1) but in analogy with the section 2.1 of the 
Matthias Blau’s lecture notes on GR, we have preferred to work with the simpler Lagrangian. 
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3 Dynamics 


We wish to construct the dynamics of the gauge field A fl and (j) such that the theory admits 
two regimes: 

1. EH regime: strong gravity reproduces the Einstein-Hilbert gravity, 

2. MOND: very weak gravity reproduces MOND and gravitational lensing attributed to 
a dark matter. 

To ensure the existence of the EH regime, we suggest that A^ and (j) are composite, 

A^ = A\ + A 2 , (3-1) 

4> = (j) 1 + </>2 • (3.2) 

We assume that the strong regime of gravity holds 

Af, = A, (3.3a) 

4> = cte, (3.3b) 

So the orbits of particles derived from (2.8) in the EH regime are identical to those derived 
from 

S[x(t)] = -^me 2 ^ J dTg lll/ x il x v , (3.4) 

which is the effective action of a particle with the effective mass of me 2 ^ in the Einstein-Hilbert 
gravity. Due to the gauge symmetry of the theory and the freedom in the mass definition, we 
choose the following representation for (3.3): 

An = 0, (3.5a) 

4> = 0. (3.5b) 


We assume that the dynamics of the metric is given by the Einstein-Hilbert action: 

S[g] = J d 4 x^-detgR. (3.6) 

In order to write the matters’ action we consider a set of particles m*, i £ {1, IV}. The world¬ 
line of these particles can be parameterized by a global choice of time in a smooth space-time 
geometry. The cumulative action of these particles, followed from (2.8), reads 


S M = dt 


' N N \ 

^2 -mi{e 2 ^gx u - 1) - ^ nmA^ 

\ i =1 i= 1 / 


(3.7) 


where the non-gravitational interactions between the particles are neglected for sake of sim¬ 
plicity. Also notice that an appropriate constant term has been added to the action in order 
to get the energy-momentum tensor of pressure-less fluid in the continuum limit (3.17). In 
the continuum approximation, (3.7) is mapped to 

SM[g^,A fl ,(j),u fl ,p\ = J dAx'J- det g (^ -^p{e 2 ^g IIIJ u IJ 'u u - 1) - A^pu^ , (3.8) 
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where p is the density and u^ is the four vector current of the matter distribution. Notice 


that due to (2.9), it holds 

e^g^u^vy = 1, (3.9) 

which should be imposed at the level of the equations of motion. 

We set the action of the gauge fields to 

S[ A i] = j d 4 x V~ det , (3.10a) 

S[A 2 ] = 16 J 9G/2 J d 4 x^d^~g£(- l ^ u Fn, (3.10b) 

where F^ u and F^ u are the gauge held strengths: 

F lu = d^ A l ~ dv A l i (3.11a) 

Ff w = d t ,Al - d v A\ , (3.11b) 


and l is a parameter, k± and k 2 are the coupling constants, and £ is a functional with the 
following asymptotic behaviors: 

r'l \ - — = j 1 ’ x > 1 
dx \ xz , 0 < x < 1 

where x > 0 is understood. We shall see that x 
represents the MOND regime of the theory. We set the action of the scalars to 

S[<f> i] = ~ 16 ^ kiG J d4x V- dets , (3.13a) 

SM = 16n l 2 Qi 2 J d A xyJ- det g C^—d^ 2 <9^ 2 ), (3.13b) 

where R\ and k 2 are the coupling constants, l and £ are introduced in the action of the gauge 
fields. 

The equations of motion of the gauge fields are derived from the variation of (3.8) and 
(3.10) : 


1V1U1NJJ regime 


represents the EH regime while x < 1 


V U F^ = IGttkiGJ^, (3.14a) 

V, (c\- l -Fl p Ff)F^ = 167 tk 2 GJ» . (3.14b) 

where = pu^ and £{z) = dC ' d A ■ The equations of the scalars follow: 

□<^>i = WnRiGp, (3.15a) 

Vj, U’(~= 16vric 2 Gp, (3.15b) 
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where (3.9) is utilized. 











In order to write the equation of the metric we first compute the energy-momentum 
tensor: 


2 2 

r J-'nu — Tmp + Ta afll' £ 
a=l a=l 

rp _ 2 <5Sm _ 2 <j> 

J-Miii' — / - j ^ ^ UfiUv i 


TA lf iu = - 

= ~ 


y/— det g Sg^ 

2 5S[Ai] 


1 


V^ddt^ 8g^ ~ 16t rmG V F ^« 


<5S[A 2 ] 


1 


\/— det g 5g^ u 1 Qttk 2 G 

2 5S[cj) i] _ 1 

° 1 ;/I/ \/ — det 5 hg /il ' 167 tkiG 

2 <55[02l 1 


/ / ;2 I 1-^211 


£{l 


)F 2 «F 2va -^X(l 2 ^) 


Tfanv — 


dpdpid^x - 2 g fll 


(3.16) 

(3.17) 

(3.18) 
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(3.19) 

(3.20) 


2 (55[0 2 ] 1 f r , n 2\9<h\\ a j. a j. 1 ^/,21<902| ^ 

= 16^iG l £ (/ ~ T^ C{1 —■ 7 


where 


|d </> 2 | 2 = -d^fad^fa, 

ra = -vr. 


(3.21) 


(3.22) 

(3.23) 


Notice that (3.9) is utilised to simplify (3.17) and write in the form of the standard energy- 
momentum tensor of pressure-less fluid: 


Tmiiv = PeffU U u fl , (3.24) 

where 

p eff = e 2 ^p. (3.25) 

The equation of motion of the metric then reads 

G — 8tt GT ^. (3.26) 

The Einstein-Hilbert regime of (3.12) simplifies the equation of motion of the gauge fields 
(3.14) and the scalars (3.15) and yields: 

V„F^ = 16tt(ki + K 2 )GJ fl , (3.27) 

□</> = 167t(ki + k 2 )Gp , (3.28) 

where — d u A^. The consistency with (3.3) demands that 


Ki T k 2 = 0 , (3.29a) 

+ k 2 = 0 . (3.29b) 

Now consider a gauge field and its held strength F /t „ and scalar (f> that solve 

= 16vrGJ^ , (3.30) 

□3> = 16nGp. (3.31) 
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After imposing appropriate boundary conditions on A\ and A 2 , and on 0 1 and 02 , one can 
writes 


A^ = , 

A% = k 2 A 9 , 


0i = , 

02 = , 


utilizing which returns: 

TAxilv + Ta 2 )iv = (al + H2 )T^j, = 0 , 
T<t>iiiv + Tfap v = (hi + k 2 )T® u = 0 , 


where 


rp.A 

1 | 

f XT' 


167 tG 1 

I r /i r VOL 

rp$ 

1 | 

{d»*d v * - 


167 tG 1 


■ v) • 


(3.32a) 

(3.32b) 


(3.32c) 

(3.32d) 


(3.33) 

(3.34) 


(3.35) 

(3.36) 


So the equation of the motion of the metric from (3.26) indeed converts to: 

G>„ = 87 tGT™ , (3.37) 

which is the Einstein-Hilbert equation for the metric. This proves that imposing (3.29) 
in the EH regime of (3.12) at the classical level reproduces the Einstein-Hilbert gravity, 
where the gauge fields and scalars affect neither the space-time geometry nor the orbits of 
particles / probes. 

Eq. (3.29) implies that one of the gauge fields abs one of the scalars have the kenitic 
energy term of ghosts. The theory, however, is equivalent to the Einstein-Hilbert theory. To 
see this define: 

A ± = ^(A] l ±A 2 ), (3.38) 

0 ± = ^(0 1 ±0 2 ), (3.39) 

Rewrite the action in term of A± and 0 . 

+ S[A 2 ] = j d^xy^detg A~d^F + ^ , (3.40) 

Sy>i] + Sy>2] = J d 4 xyj- det g 0 _ D 0 + , (3.41) 

The exact Feynman path integral (up to a normalisation factor), then, is simplified to: 

Z = J DgDA + DA-Dct) + D(j)-e~^ (3.42) 

= [ DgDA + D(t> + 8{d il F+ u )8{U^ + )e-^ s[9]+SM) , 


(3.43) 








where the integrations over A and 0 are performed. The simplified path integral gives the 
exact equations of A + and 0 + at the level of quantum held theory: 


which are solved by 


d' J F+ u = 0, (3.44) 

□0+ = 0, (3.45) 

A+ = 0, (3.46) 

0 + = 0, (3.47) 


inserting which into (3.42) reproduces the exact Feynman path integral of the Einstein-Hilbert 
action. This proves the EH limit of the theory coincides to the Einstein-Hilbert theory at the 
level of quantum held theory too. 


4 Static Solutions 


Let a static mass distribution be considered in an asymptotically hat space-time. The metric 
in the not-strong regime of the Einstein-Hilbert part of the action reads: 

ds 2 = —(1 + 2 cf>]y)dt 2 + (1 — 2 <f>N)dx 2 , (4-1) 


where holds 

and p is the density, and 2 <f>N -C 1. 


V 2 0jv = A-nGp, (4.2) 

The gauge helds and scalars can be represented by: 


A llx = (A^x), 0,0,0), (4.3) 

A 2 » = (A 2 (x), 0,0,0), (4.4) 

0i = 0i(x), (4.5) 

02 = fa(x ). (4.6) 


The four vector velocity held of the matter distribution reads: 


= (u(x), 0,0, 0), 


(4.7) 


Notice that (3.9) then implies 


u(x) 


e -0O) 
i + 4>n 


where 0(x) = <f>i(x) + 0 2 ( 3 ;). Utilizing this metric in (3.15), and using the identity of 


(4.8) 


□0 = 1 d^-det gg^d,,#) , 

V — det g 


(4.9) 


0 1 = 4&107V , 

02 = 4fc 2 07V , 
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and recalling that 0 jv 1, returns 


(4.10) 

(4.11) 






The equations of the gauge fields in the Einstein-Hilbert regime gives: 


A\ = 4:k\(j) N , (4.12) 

A'2 = 4:k 2 (f > N , (4.13) 


where the boundary conditions are imposed such that the fields and their first derivatives 
vanish for vanishing total mass. 

Although the net contribution of the gauge fields and scalars to the energy momentum 
tensor does not vanish in the MOND regime, it is negligible assuming. The metric in the 
MOND regime, therefore, coincides to (4.1) at the leading order approximation and can be 
approximated to the flat Minkowski metric. This simplifies the equation of motion in the 
MOND regime -(3.14b) and (3.15b) for static solutions to: 


v‘(®4v*) 

V( ^v,A 2 ) 


16irk 2 Gp , 
16irK 2 Gp , 


Recall the AQUAL equation in its MOND regime [2]: 

V 4 =4v r Gp. 

So we can write: 

(j) 2 (r) = 2sign(ic 2 )^0 m , 

A 2 {r) = 2sign(/i 2 )\/M 0m , 

where: 

c 2 

a* = —7=- . 

V21 


(4.14) 

(4.15) 


(4-16) 


(4.17a) 

(4.17b) 


(4-18) 


4.1 Orbit of particles 

Eq. (2.8) that governs orbits of massive particle for static mass distribution simplifies to: 


5[x(r)] = J drL (4.19) 

L = 4-—— (—(1 + 2</>jv)i 2 + (1 — 2 <^>tv) |ar| 2 ) — mA^i . (4.20) 

For a slow moving massive particle, one can choose 

t = t, (4.21) 

|f| < 1, (4.22) 

that simplifies (4.19) to 

S = f dt (^m\x\ 2 - m(A 0 + cj) + c^n)) + 0(4> 2 , <f> 2 N ) . (4.23) 
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So the effective weak gravitational potential for a slow moving probe follows 

fieff 1 = + 4> + 4>n ■ (4.24) 

Now let a fast moving particle be considered. We can choose the following coordinate 


t = t , 

(4.25a) 

x 3 = t + x 3 , 

(4.25b) 

Xl = Xl, 

(4.25c) 

X2 = X 2 , 

(4.25d) 

with the world-line parametrization of (4.21) such that 


\x\ , |s 2 |, |£3| < 1 • 

(4.26) 


Inserting (4.25) into (2.8) gives 

S = j dt (^m\x\ 2 - m(A 0 + 2 <j>)) + 0(</> 2 , 4>n) ■ (4.27) 


which means that the effective weak gravitational potential for a fast moving probe including 
a light signal reads 


0 eff t — Tlo + 2 (/>jv ■ 


(4.28) 


Notice that as expected, (f> does not affect the light’s trajectory. The Newtonian regime due 
to Aq = (j) = 0 holds 


,fast 

p eff 

■k,slow 


= 7 + 1 = 2 , 


(4.29) 


Veff 

which means that the post Newtonian parameter of 7 = 1. Since the MOND regime also 
holds 7=1 [5] then we conclude from (4.24) and (4.28) that the MOND regime should hold 


(/) = 



(4.30) 


which means that the scalars decrease the value of effective gravitational potential to 50%. 
We notice that (4.30) can not be hold in the interpolating regime: the regime that connects 
the Newtonian regime to the MONDian regime. We discuss the MOND regime in the next 
section. Einstein-Hilbert regime holds (4.30). In the MOND regime, after approximating Aq 
and 0 respectively to A 2 and cj) 2 , (4.30) demands 


k 2 = -^« 2 , (4.31) 

which beside (3.29) allows to express all the coupling constants in term of a constant: 

k j = —k , (4.32) 

k 2 = +^, (4.33) 

hi = (4.34) 

k 2 = (4.35) 
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Far from the mass distribution holds 


d r (j)N = + (4.36a) 

\JGMa * „ ^ 1 ^ 

<9 r </> m = - r + 0(-z), (4.36b) 

iy* j • — 

where M represents the total mass. Inserting (4.36) into (4.10), (4.12) and (4.17), and using 
the results in (4.24) gives 

dr€ l f7 = sign (fc) yggjK + 0(1) . (4.37) 

We notice that requiring the gravitational attraction at large r demands that 

k>0. (4.38) 


In order to reproduce the flat rotational velocity curves of spiral galaxies, the effective poten¬ 
tial should coincide to [ 2 ]: 


dAiff' = 4^± + o(±), 

a Q = (1 ± 0 . 2 ) x 10" 10? 5 . 

s z 

The consistency between (4.18), (4.37) and (4.39) then demand: 

l _ c 2 
K y/2 a Q ' 


(4.39a) 

(4.39b) 


(4.40) 


5 Weak Gravitational Lensing 

Since l is (4.40) sufficiently large, we can approximate (4.1) to the Minkowsky metric and 
simplify (3.15b) and (3.14b) around a spherical static mass distribution to: 

12 CM 

£'(—\d r A 2 \ 2 )d r A 2 = 4k^- , (5.1a) 

/ 2 CM 

= -K-^2- , (5-lb) 

where M is the total mass. To find the exact value of the fields we consider the following 
family of the interpolating functions: 

C>(x 2 ) = f, n (x) = - — —j- , (5.2) 

(l+X n )n 

where x > 0. We can write 

£ n (x) = 2^1 ( — ,1 ± ~x^) + Cl , (5.3) 

3 n n n 

where 2 F\ is the Gauss’s hypergeometric function and ci is the constant of integration. Notice 
that 

lim C n (x) = ci, (5.4) 

X-»0 
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and ci is tantamount to adding a cosmological constant. For the time being we set 


The first few C n {x) follows 



ci = 0 . 

(5.5) 

£i(x) 

= —2y/x + x + 2 ln(l + v 7 ®), 

(5.6) 

£3 (x) 
2 

= y 7 x 3 / 4 + l(x 3 / 4 — 3) + 3 , 

(5.7) 

£ 2 ( 2 ;) 

= y 7 x(l + x) — ArcSinh(v 7 x), 

(5.8) 

£ 3 ( 2 ;) 

= -! + (! + x§)i , 

(5.9) 


The literature knows the cases of n = 1 and n = 2 respectively as the simple and the standard 
choices. Define new variables 


l 


Xn. = 


x<t> = 


V2c 


d r A 2 , 


V2c 2 


d r <t> 2 


and insert them into (5.1): 


b <t> 


x 2 

(i + ®S)" 


where (4.40) is ultilized and 


__ 1 GM 


(5.10) is solved by: 


nr I ,Kr,„ 

H = 2 y { 2 = (—) + I/ + <“) 


i-i nr ( n ,nr ,„ I „ . nr . n 

*« = 2 8U-)" + t/4 + 8»(-)2» 


which utilizing (5.10) gives 


Eq. (4.10) and (4.12) give: 


,/ _ «o 

</>2 — 2k X ^’ 

A' - —T 

2 - 2 k X q . 


J — 

<^1 = Ka 0 - 2 

^2 

A) = -4«:a 0 ^ 


(5.10a) 

(5.10b) 


,7.2 

- 2 k 2 

X r z 

(5.11) 

r 2 

- = 8k 2 ^ , 
r r z 

(5.12) 


(5.13) 

(5.14) 

(5.15) 

(5.16) 

(5.17) 

(5.18) 

(5.19) 
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Jn- 2 in~ 2 



f."2 in ~2 



in- 2 in~ 2 



Figure 1: The 7 for various values of k and different choices of the interpolating function. 


Now let 7 be defined by 


r) ch^ ast 
°r9 e ff 

f) AsloW 
u r cp e ff 


t 4 q + 2 (f/ N 

A' 0 + 4>' + 4>' n 


(5.20) 


in the 1. h. s. of which (4.24) and (4.28) are utilized and the prime superscript represents 
differentiation with respect to r. The Einstein-Hilbert gravity requires that 7=1, and in or 
study it is easier to work with 7 rather than 7 . The exact form of 7 for any choice of /r(n) in 
(5.2) follows: 
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'Slow 


0Fasl 

(IQ ClQ 



do do 



a 0 d 0 



Figure 2: The effective gravitational acceleration felt by a slow (left) and fast (right) moving 
probe for various values of k and different choices of the interpolating function. 


-2 1+ s(-l + k) + & ( ^4 + 2" (f ) 2n + 2-/ 2 (f ) n ) n 

— = 1 +-^-^- T . (5.21) 

7 ” 2 2+ ^(-1 + 2k)- ^ ^4 + 8 -(f) 2n + 2 3-/2 (f)^" 


Fig. 1 depicts 7 n for various values of n and k. As depicted in the fig. 2 , the sharp increase 
in 7 is due to the sharp decrease of 4>ffJ■ We see that requiring 7 to be close to 2 sets 
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Figure 3: The colored regimes shows the values of k and n that are consistent with the weak 
gravitational lensing at one, two and three sigma confidence levels. 


non-trivial constraints on k and n. Ref. [5] reports that 

I 7 — 1 | 

1 ' 1 < 0.02, (5.22) 

at one sigma level. This means that 7 is constant at the leading order approximation. There¬ 
fore, it equivalently holds 


17 n 2 

2 


< 0 . 02 . 


We notice that the maximum value of 7 holds 7 ' n (r) = 0. So 


7 n(r) = 2 +CL x 0.04, 
7n( r ) = 0 


(5.23) 


(5.24) 

(5.25) 


can numerically be solved in order to find the maximum value of k that holds (5.23) everywhere 
with the confidence level of CL. Fig. 3 depicts the values of k and n holding (5.23) everywhere 
at one, two and three sigma confidence levels. For these allowed values, no sharp decrease 
exists in the effective potentials. 

The 7 parameters is measured in strong gravity too. One of the best bound on its strong 
gravity value reads [ 8 , 9]: 

7 - 1 = (0.8 ± 1.2) x 10“ 4 , (5.26) 

Fig. 4 depicts the value of 7 n at the surface of the Sun for various values of k and n. We 
observe that they are in agreement with (5.26). 


6 Lunar System Constraints 

The effective GM is defined by 

GM eff (r) = r 2 d r (j) s ^f. 
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( 6 . 1 ) 












7,1 -2 



Figure 4: 7 n — 2 at the surface of the Sun for various values of n and k. 


The accurate value of the mass ratio of the Sun/(Earth+Moon) from the Lunar Laser Ranging 
combined with the Solar GM and the lunar GM from lunar orbiting spacecrafts [6] yields 
the effective gravitational mass of the Earth in an Earth-centered reference frame with the 
precision of one part in 10 8 : 

GM ELR (r LD ) = 398600.443 ± 0.004—(6.2) 
where r LD represents the Lunar distance: the average distance between the Earth and Moon: 

t ld = 384,400/cm , (6.3) 


The Earth GM has also been measured by various artificial Earth satellites, including the 
accurate tracking of the LAGEOS satellites orbiting the Earth in nearly circular orbits with 
semimajor axes about twice the radius of the Earth [7]: 

GM RAGEOS (2r (B ) = 398600.4419 ± 0.0002-2- (6.4) 

where r® stands for the radius of the Earth: 

r® = 6, 371km . (6.5) 


The consistency with the Lunar system, thus, requires 


GM e J f (r LD ) - GM e J f { 2r®) 
GM® 


< 10 " 8 . 


( 6 . 6 ) 


We notice that the distance between the Moon and the earth are measured by sending light 
signals. Deviation from the Einstein-Hilbert geometry affects the time that a signal travel 
from the Earth to the moon and back, cj) does not contribute to null geodesics. Our theory 
adds AT to the time of travel for the signal: 


2 f T LD 

AT = — / dr At 

C " 4r e 




24>n) , 


(6.7) 


- 16 - 








wherein we have approximated dt = ^dr since the space-time geometry around the Earth is 
almost flat. We can taylor expands the fields in the Newtonian regime of the Lunar system,, 
in the regime that holds ( —) 2 l 2 : 


9r€ l f7= „2 

8 d>f ast — 

Ur<P e ff ~ 


GM 


r 

2 GM 


1 _ 2 2n — 4 n nr 2 n 
2 3n rGr 


+ 


"eff 


n r 


+ 


In order to have a consistent Taylor series for the Lunar system, we demand that 

(-*^)- 2?l < 10~ p ,-> 
v8 r LD 

K\(n,p) = 8\/2 X 10^ -4 , 

Ki(n,p) < k, 

where p > 0 and the larger p is the better the Taylor series is. Eq. (6.9) leads to 


cAT = 


}1—3n^.l— 2n 


^ ^LD ^2 n ^ w 


2 r, 


© \ 2n 


GM a 


n 2 (2n - 1) V r ~e ' r ® ' J <? 

wherein ~ 20tv is assumed to vanish at r = 0. Notice that 

^ = 4xl0- 4 , 

r® 

= 1.2 x 10 -5 , 


GM($ 9 o 

—^ = - x 10" 3 m. 
c z 1 


( 6 . 8 ) 

(6.9) 

(6.10a) 

(6.10b) 

(6.10c) 

( 6 . 11 ) 

( 6 . 12 ) 

(6.13) 

(6.14) 


The LLR measurements identify the center of the earth with the precision of few centimeters. 
We require 

cAT < 1cm , (6.15) 

which induces the following lower bound on k: 


K 2 {n = 


9 x 2 r> 


n 2 (2n — 1) 


i 

2n —1 


10 


-4-7 


K 2 (n) < K . 


(6.16a) 

(6.16b) 


This lower bound allows us to neglects the corrections of the theory to the light speed 
travel time between the Moon and Earth. Notice that when (6.15) is not met, one should 
repeat the calculation of the Earth GM mass from raw data within our theory. This is 
clearly computationally extensive. We have escaped from this task by considering the set of 
parameter that holds (6.15). We further observe that the lower bound in (6.10a) for all p > 0 
is larger than the lower bound in (6.15) for n > 1. The consistency with the Lunar system, 
then, requires (6.6). Utilizing (6.9) then simplifies (6.6) to 


K 3 (n) = 

K 3 (n) < K 


(2 2n - 4)2^ 


n 


10 


— 4+7 


2 n = 1 does not resolve the Pioneer anomaly. 


(6.17a) 

(6.17b) 
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(b) Three sigmas. 

Figure 5: The values of n and n which are consistent with the Lunar system data for 
p = 1, 2, • • • ,6 and the weak gravitational lensing at one and three sigmas . The smaller the 
regime is the larger p is. 


Notice that the lunar system provides the lower bounds on k while the gravitational lensing 
provide an upper bond. The fig. 5a shows the values of n and k that are consistent with (6.6) 
and (6.10a) and one sigma level constraint in fig. 3. The allowed regimes are disconnected. 
The arbitrary large values of n are allowed for k < 0.02. The fig. 5b depicts the values of n 
and k that are consistent with (6.6) and (6.10a) and the three sigma level constraint of fig. 
3. 

7 On the consistency of the static spherical solutions 

At the end of section 3, we have shown that the energy momentum tensor of the gauge fields 
and scalar vanish in the limit of the Einstein-Hilbert gravity. In this section we would like to 
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compute the correction to the energy momentum tensor as we move away from the Einstein- 
Hilbert regime. In so doing, we notice that the following taylor expansions exist around the 
Einstein-Hilbert regime: 


/ \ i 1 — n 2 1 + n _2 n 4 Id - 3 n + 2 n _o„ g 

Hn{x) = l--x e 2 + — 2 e 4 - —7 - x 6n e b + • 

n 2 n z or?. 


£„(x 2 ) = x ( 1 — 


(2 — n)n 


—n 2 1 1 + u _ 2 n 4 1 + 3n + 2 n _o n g 

x e 2 + —- -^x 2 e 4 - 7 -- -x 6n e b + 


2(1 — n)n 2 


3n 3 (2 - 3n) 


(7.1) 

(7.2) 


wherein the auxiliary parameter of e is turned on to systematically track the perturbation 
series. At the end of the perturbative computation we set e = 1. We also get: 


167tG y (T4 |0 o + Ta 2 oo) — 77o ( — 


c 4 , 

( kf 

kl' 2 1 

\ r 

c 4 , 

( kf 

kl' 2 1 

\ r 


4-2n / 2 7 “ 3n 


4—2n 


n(n — 2) 
2 5 ~ 


e 2 - 


2 5 - 6 n( 3n _ ^ /^\ 


— 2n 


167 i-G(r*oo + T fc oo) -^ 2 ( 7 -) l-7irT7 e2 + d?7 : VT7 I ( : 7j 


,, c (7.3) 

n 2 (n — 1 ) \ r ) I 

2 3 _ 2 n (3n — 1) / kf\~ 2n 4 \ 


(7.4) 


So 


16ir0r » = «dv) 


Let it be dehned: 


4—2n 


2 5 - 3 n (4 n - 4) 2 (4 2n - 4)(3n - 1) f kf\ 

(n — 2)n € 8 2n_1 (ro — l)n 2 \ r J 


—2n 


M eff( r ) = 4 / ^00 ( r 0 4:irr ,2 dr' 
c Jo 


(7.5) 

(7.6) 


M e ff( r ) is the effective mass that Too effectively induces in a sphere of radius r around the 
origin. We notice that only for 

\<n, (7.7) 

the effective mass is bounded. We further notice that for 


- < n < 1, 


and 


2 < n , 

the effective mass around the Einstein-Hilbert regime is negative while it is positive for 


(7.8) 

(7.9) 


1 < n < 2 . 


In the MOND regime we have the following perturbative series: 

Ah i(x) = Vx( 1 - -x?e 2 ), 
n 

. 2 3 3 — 2\ 

£»(x) = ^ 2 l - 7 X— 

3 n(3 + n) 


(7.10) 

(7.11) 

(7.12) 
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The perturbation results to 


167tGToo = 


48c 4 f kr\ 4 32\/2c 4 [hr 


- 


kl 2 V 




3 kl 2 V r 


1 + 3x22 


_ 1 (2 n+1 - 1) fkr 


n(3 + n) \ r 


(7.13) 


We observe that for all values of n > 0, the tree-level Too is negative in the MOND regime 
for: 


9 r 

--K < - . 


2y/2 r 


(7.14) 


Now let it be defined 


M e ff(r) = ~2 Too(r') Airr^dr 1 , 


(7.15) 


Meff(r) is the effective mass that Too effectively induces between a sphere of a radius r and 
a sphere of radius L. For a sufficient large value of L we get: 


v^2 c 2 ^ 2 3 L_ 

3 GP [ ’ [ kr h 


Meff = ~ ^ V3 


(7.16) 


For | M e ff | at the order of M then the perturbative solution around the flat geometry is not 
valid. This means that there is a minimum value for L for which the perturbation is valid. 
We can find its order by demanding M e ff = —M: 


M 0 = = 1.21 x 10 54 %, 

a 0 G 


L 

2 K 





(7.17) 

(7.18) 


The hnite value for L means that the MONDian £ behavior does not extend to infinity. At 
order L, the gravitational held must declines faster than A For the Sun, L reads 


=3.3xl0 71658595155 m. 

K 


(7.19) 


For a galaxy of mass M g = 10 12 M@, we get 

^ = 4.5 x 10 720216 ?n . (7.20) 

Hi 

For a clauster galaxy of mass M c = 10 16 Mq, it yields 

— = 4.9 x 10 7222 m. (7.21) 

Hi 

Due to the size of the Universe we can not observe the deviation from - in the deep MONDian 
regime. 
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Figure 6: The values of n and n that are consistent with 3<r weak gravitational lensing, 
Lunar system and (8.8). 

8 Cosmological Constant 

Assuming that the constant of integration is not vanishing in (5.3) leads to the following term 
in the action: 



( 8 . 1 ) 


which is equivalent to adding a cosmological term to the action: 



( 8 . 2 ) 


The cosmological constant in term of h and 12a reads 


A = 4 Hfo A . 

c z 


(8.3) 


Table IV of [10] reports values of 12 a and Hq due to +SN la observation as : 


12 a = 0.725; 
h = 0.599 


+ 0.039 
'- 0.044 5 


(8.4) 

(8.5) 


.+ 0.090 
- 0.062 ) 
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Figure 7: 7 n (r) for n = 1 and k = g and g. 


where Hq = 100 h km s 1 MPc 1 = 3.24 x 10 m hs 1 . Using these values, one finds 


—18 u D—i 


Ameasured = 1.206±2‘g^ X 10~ 52 —- 


meters 


2 ■ 


Equating (8.6) to (8.2) identifies C\: 

K C A m e aS ur e d 


Cl = 


6 al 


= (18.70 ± 8.08) At 3 


We demand ci and k to be at order one. This leads to 

0.15 < k < 0.70. 


( 8 . 6 ) 


(8.7) 


( 8 . 8 ) 


Fig 6 the allowed values of n and k consistent with (8.8) and fig. 5b. This demands that 
n = 1. We prefer a value for n equals to the inverse of a natural number: 


and 


K 

n 

1 

ci 


1 

5 ’ 

1, 

5 • ••11 


(8.9a) 

(8.9b) 

(8.9c) 


k = 1 , ( 8 .10a) 

6 

n = 1 , ( 8 .10b) 

— = 8---20 ( 8 . 10 c) 

ci 


reproduce the observed value of the cosmological constant, generates MOND and the weak 
gravitational lensing attributed to dark matter at three sigma, and are consistent with the 
lunar system data. Fig. 7 presents 7 for (8.9). We notice that the average of A 7 in E E (0,10) 
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for re = | reads: 
5 


7-2 

< >= 


r-10f 


f (lOf ) 3 A, 


dr47rr 2 


7-2 


= 0.014 


while for re = J it reads 


7-2 

< >= 


f (I0f) 3 7 0 


10? o' — 2 

drAnr 2 - = 0.017. 


( 8 . 11 ) 


( 8 . 12 ) 


This means that these values of re are in fact in agreement with [5] that reports the average of 
7. We observe that only for 1 £ (0.15,6), Ay = exceeds 0.02. In particular the average 
of Ay in 1 £ (^,2) exceeds well beyond 0.02: 

fc = i -» < Ay >= 0.05, (8.13) 


5 

k =-> < Ay >= 0.04. 

6 


(8.14) 


This means that we should find the average of < Ay > for 1 £ (f, 2) rather than its average 
over whole of the space. This is a prediction of the theory. This prediction can be attested 
by refining the current data or specifically designed future observations. 

9 Cosmological Solutions 

We accept that cosmology is described by a homogeneous isotropic solution in the sense that 
the the large matter inhomogeneities in small scales do not produce significant effects in large 
scales [11]. 

The field content of a homogeneous isotropic time dependent solution of the theory is 
represented by 

ds 2 = — dt 2 + a(t) 2 d\x\ 2 , (9.1) 

fa = fa(t), (9-2) 

fa = fa(t), (9-3) 

4 = A\ = (0,0, 0,0). (9.4) 

Notice that the gauge fields vanish due to symmetries. We represent the local matter density 
by p = p(t), and its four velocity by: 

(9.5) 

(3.9) implies that 


U/i = K 0,0,0) , 
e 2 ^Ug = 1 —> uq = e~ 


(9.6) 

(9.7) 


So the matter current reads: 

J tl = pe~*( TO, 0,0). 

We observe that a vanishing gauge field does not satisfy (3.14). 3 We need to change our 
theory to become consistent with cosmology. To do so define: 

f d 3 x\/— det g 




f d?xy/— det g 


(9.8) 


’This paradox also exists in ordinary QED: The electric field of a uniformly charged space vanish due to 
symmetries but the vanishing field is not consistent with the gauss law. 
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Next alter the matter action (3.8) to: 


Sm = J d 4 xy/~ det g(-)^pe 2(,, g^u ,J 'u u - A^J 11 - J»)) (9.9) 

The equations of motion of the gauge-helds then follow 

V v F{' 11 = 16ttki G(J^ - J 11 ), (9.10a) 

V, = 167tk 2 C(J^ - J ' 1 ). (9.10b) 

Since any local mass distribution in a non-compact universe holds = 0 then the solutions of 
the previous sections are left intact. The homogeneous isotropic cosmological solution (9.7), 
however, holds: 

= 0. (9.11) 

The equations of motion of the gauge-helds (9.10), therefore, are in agreement with a vanishing 
gauge fields (9.4). 

In order to address the dynamics of 02, we need to identify C{x) for negative value of x 
because 


- <^02(9^02 = -02 (t) 2 < 0 

(9.12) 

We consider the simplest choice: 


C{x) = x + c\ Vx < 0 . 

(9.13) 


Notice that ci should be present in order to have C in C 1 . The equations of motion of the 
scalars then read: 


= 167 tkiG/9 , (9.14a) 

□02 = 16tvK 2 Gp , (9.14b) 

Due to the form of these equations, (3.29) and choosing appropriate boundary conditions we 
get 

0 = 01 + 02 = 0 . (9.15) 

This means that the scalars as well as the gauge fields do not affects the trajectories of the 
particles. The net contribution of the energy momentum tensors of the scalar and gauge fields 
vanish too: 


+ T4>2hu — 0 , (9.16) 

Ta^u + Ta^u = 0 . ( 9 - 17 ) 

Notice that the ci in the definition of L generates the cosmological constant. The equation 
of motion of the metric then simplifies to: 

+ Ag jiv = 8ttGTmi_il' 1 (9.18) 

where A is given in (8.2). So the cosmological solution of the theory coincides to that of 
the Einstein-Hilbert action. The late-time behavior of the theory is in agreement with the 
observation. 
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10 Conclusions and discussions 


This research has been started by asking how precise the space-time geometry can be ap¬ 
proximated by a pseudo Riemann geometry. Can quantum gravity or dark matter/energy be 
naturally resolved in geometries beyond Riemann? Can the degrees of freedom for dark fields 
be assigned to a deviation from Riemann geometry instead of being added to the table of the 
ordinary elementary particles? Can the MOND theory [1] be naturally realised in a deviation 
from the Riemann geometry? 

In order to answer some of these questions, we have considered a deviation from the 
Riemann geometry toward the Randers’ geometry in the section 2. This deviation includes 
two gauge fields and two scalars in addition to the metric. It possesses four coupling constants, 
a length scale and a functional. 

In the section 3, we have fixed the dynamics of the gauge fields such that the strong 
field limit of the theory coincides to that of the Einstein-Hilbert theory. This has reduced 
the number of the free coupling constants to two. We have further fixed the dynamics of the 
theory such that its weak field limit reproduces the AQUAL [2] and MOND theories. This 
has mapped the length parameter of the theory to the critical acceleration of the MOND, and 
has fixed the asymptotic behaviours of the functional of theory. 

In section 4, we have presented the static solutions in the MOND and Newtonian regimes 
of the theory. We have studied the light trajectory around static solutions. We have calculated 
the 7 parameter in the MOND regime of theory. We have required the 7 parameter in the 
MOND regime to be one. This requirement has reduced the number of coupling constants to 
one. 

We have considered a family for the functional of the theory in the section 5. This 
family can be expressed in terms of the Gauss’s hypergeometric function and a constant 
of integration. We have assumed that the constant of the integration is sufficiently small 
and does not affect the space-time geometry at the considered vicinity of the spherical mass 
distribution. We have ignored the negligible contribution of the gauge fields and scalars in 
the MOND regime to the metric by their contributions to the energy momentum tensor. In 
doing so, we have found the exact spherical solutions of the theory. We have calculated the 
7 parameter in the whole of the space. We have reported that requiring the 7 parameter to 
be consistent with the constraint reported in [5] sets a strong upper bound on the coupling 
constant of the theory as depicted in the fig. 3. We have observed that this bound is in 
agreement with the value of 7 measured at the surface of the Sun as depicted in fig. 4. 

We have studied the Lunar System constraints on the theory in the section 6 . We have 
used the Laser Lunar Ranging measurements [ 6 ] and LAGEOS [7] data. We have assumed 
that the space-time geometry around the earth provides a perturbation to that predicted by 
the Einstein-Hilbert gravity. We have sufficed to the inclusion of the leading perturbation and 
assumed that this perturbation is smaller by the factor of 10 -p from the Newtonian gravity. 
We have identified the regime of the parameters wherein the theory does not affect the time 
of the travel of the signals from the Earth to the Moon with the precision reported in [ 6 ]. 
We have shown that these considerations and [ 6 , 7] provide a lower bound on the coupling 
constant of the theory as depicted in the fig. 5a and 5b. 

We have studied the consistency of our approach in 7. We have shown that as we move 
toward the asymptotic infinity, the neglected contribution of the energy-momentum tensor of 
the gauge and scalar fields accumulates. For any given mass, we have calculated the length 
scale wherein the energy-momentum tensor of the gauge and scalar fields must be considered, 
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the result has been expressed in (7.18). We have reported that for any given astronomical 
mass, this length scale is far larger than the size of the observed Universe. Notice that the 
kinetic term of one of the gauge fields and one of the scalars have negative sign and they 
might cause a new difficulty in the quantisation of the geometry that here is composed of the 
metric, two scalars (non-minimally coupled to matter) and two gauge fields. We have shown 
that the strong field limit of theory coincides to the Einstein-Hilbert gravity also at the level 
of quantum field theory. We have studied the classical solutions (IR limit) and we have shown 
that the theory is classically consistent in its IR limit. The theory is not trivial. Addressing 
the quantum stability of the theory in the interpolating IR-UV regime requires addressing 
quantisation of the geometry and is not within the scope of this work. 

We have mapped the constant of the integration in the definition of the functional of 
the theory to the observed cosmological constant in the section of 8 . This map is natural in 
the sense that it needs a tuning at order one. So we have shown that the theory naturally 
accommodates the cosmological constant. We further have shown that the accommodation 
identifies the only free coupling constant of the theory. We have reported that the identified 
theory predicts that the 7 parameter in the interpolating regime deviates from 1 by few 
percent albeit only in a specific section of the interpolating regime. This prediction can be 
attested by refining the current data or specified future observations. 

We have studied the cosmology of the theory in section 9. We have slightly changed the 
coupling of the gauge fields to the matter current in order to have a homogeneous isotropic 
time-dependent solution. This change leaves intact the static solutions of the previous sec¬ 
tions. We have reported that the isotropic homogeneous solution of the theory coincides to 
that of the Einstein-Hilbert theory with a cosmological constant. So the late time cosmology 
of the theory is indeed in accord with observation. We have left the study of the early universe 
in particular that of the fluctuations of CMB in the theory to a future study. 

In summery, we have presented a deviation of the physical length from the Riemann 
geometry toward the Randers’ geometry. We have constructed a consistent second-order rel¬ 
ativistic theory of gravity that dynamically reduces to the Einstein-Hilbert theory for strong 
and Newtonian gravity while its weak gravitational regime reproduces MOND and the grav¬ 
itational lensing attributed to the dark matter halo. It also naturally accommodates the 
observed value of the cosmological constant. We have reported that it predicts a small de¬ 
viation for the post Newtonian parameter 7 from 1 in a specific regime that interpolates 
the Newtonian regime to the MOND regime. The deviation is consistent with the reported 
observations but can possibly be detected by fine-tuned refinements of the current or future 
data. We, however, do not jump into the conclusion that all the reported aspects of dark 
matter, including the acoustic peaks in CMB or the Bullet cluster, can be reproduced by this 
model. Since the model is mathematically simple and phenomenologically adaptable, we do 
not refute the possibility that further aspects of dark matter can be accommodated in de¬ 
grees of freedom of the space-time geometry that encodes some deviations from the Riemann 
geometry. Future works will investigate this possibility. 
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